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Abstract 


The  problem  of  direct  Interpretation  of  earth  resistivity 
data  was  attacked  by  Ness  (4)  using  a  method  of  least  squares. 

The  same  problem  is  considered  here,  using  linear  programming 
methods  for  its  solution.  These  methods  were  tried  because 
they  allow  data  fitting  in  other  than  a  least  square  sense. 

Two  sets  of  model  data  were  analyzed,  and  the  results 
are  comparable  with  those  obtained  by  Ness,  both  in  quality 
of  fit  and  in  computation  time  required.  It  appears  that  both 
methods  reduce  errors  to  within  the  range  of  the  accuracy  of 
the  measurements  and  of  the  approximations. 

The  linear  programming  methods  possess  a  decided  advan¬ 
tage  in  their  flexibility,  which  allow  subjective  control 
where  desired,  and  which  would  also  allow  the  use  of  informa¬ 
tion  from  other  sources  with  little  modification  of  the  programs. 
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I .  Introduction 


The  direct  interpretation  of  resistivity  data  requires 
the  solution  of  an  inverse  boundary  value  problem.  The  main 
work  in  this  connection  was  done  by  Vozoff  [6,7]  and  Ness  M. 
In  [4]  one  finds  a  considerable  historical  account  of  resis¬ 
tivity  measurements,  and  of  attempts  at  direct  and  indirect 
interpretation. 

The  direct  interpretation  problem  has  been  formulated 
by  a  first  approximation,  as  a  set  of  linear  equations  and 
adapted  for  computer  solution  M  [6,7].  The  method  of  least 
squares  was  used  in  fitting  the  observed  data.  It  was  indi¬ 
cated  that  this  method  gave  results  which  may  be  classed  as 
semi quantitative .  Detailed  discussion  on  this  aspect  of  the 
problem  is  found  in  [4,6, 7D  * 

Linear  programming  is,  comparatively,  a  new  mathematical 
method  for  solving  problems  concerned  with  decision-making. 

The  method  of  solution  is  basically  iterative,  and  thus  ideally 
suited  for  electronic  computers.  Further,  the  flexibility  of 
its  methods  give  a  powerful  control  over  the  solution  itself. 

A  large  variety  of  problems  are  amenable  to  analysis  by 
linear  programming.  It  can  be  used,  for  example,  to  minimize 
the  sum  of  absolute  value  of  errors,  or  to  emphasize  or  ignore 
the  effect  of  a  number  of  variables. 
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In  the  following,  an  attempt  is  made 
of  linear  programming  in  interpreting  the 
directly. 


to  use  the  method 
resistivity  data 
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II.  Formulation  of  the  Problem 

The  problem  of  interpreting  the  apparent  resistivity  data, 
assuming  general  inhomogeneity  is  fundamentally  based  on  the 
following  hypothesis,  originally  put  forward  by  Vozoff 

"The  region  of  interest  can  be  considered  as  being  made 
up  of  small  homogeneous  blocks  of  given  geometry  but  unknown 
conductivity. " 

This  leads  to  a  three  or  two  dimensional  array  of  regular 
blocks  of  known  geometry.  The  simplest  shape  to  choose 
would  be  that  of  a  cube.  Accordingly  in  the  following  we 
shall  assume  the  subsurface  region  of  interest  to  be  made  up 
of  an  array  of  small  homogeneous  cubes  having  unknown  con¬ 
ductivities  3  i  —  lj  2-; . il  ,  where  n  is  the  total 

number  of  cubes  into  which  the  subsurface  region  of  interest 
has  been  divided.  This  matrix  of  n  cubes  is  assumed  to  be 
embedded  in  an  otherwise  semi-infinite  homogeneous  medium 
characterized  by  ,  the  background  conductivity. 

Suppose  is  an  apparent  resistivity  reading,  at  a 

particular  location  on  the  surface  of  the  earth,  or  on  the 
model  in  the  laboratory.  We  assume  the  combined  effect  of 
the  n  cubes  together  with  the  primary  effect  to  be  equal  to  ra. 
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Thus  symbolically  we  have 


n 


+  F  (<4,o;,.. .  ok,  ^  %z) 


where  F  is  some  general  function  of  the  arguments  ^  ) 


1  =  I j  2/j  •  •  •  ■  n  and  x,  y,  z. 

The  approximation  is  such  that  it  linearizes  the  system 
and  expresses  F  as 


where 


=  geometric  factor  associated  with 
X  fit  cube  and  the  relative  positions 
of  source  and  receiver. 


We  shall  state  here  the  final  results  of  the  first 
approximation,  referring  for  actual  details  to  the  references 
cited. 

The  linearized  system  of  equations  for  pole-pole  configura¬ 
tion  (see  Fig.  l)  i.e.  for  the  potential  at  a  point  P  on  the 
surface  due  to  a  unit  current  source  at  another  point  Q  on 
the  surface  is 
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GEOMETRICAL  PICTURE  FOR  A  SINGLE  CUBE 
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We  shall  use 


(a)  Index  i  to  refer  to  the  cubes  into  which  the  subsur¬ 
face  region  has  been  divided,  (i  =  1,  2y  ....n) 

(b)  index  j  to  refer  to  receiver  point  P.  (in  numerical 

calculation  values  j  =  -5y  -4*  -3*  . . 3,  4„  5 

were  given). 


(c)  index  k  to  refer  to  source  point  Q.  (in  numerical 
calculations  values  k  -  -5*  -4,  ....... 4„  5  were 

given) . 


An  apparent  resistivity  reading  is  characterized  by  a 
combination  of  j  and  k  index.  We  shall  denote  the  combination 
jk  by  m.  Thus  m  refers  to  apparent  resistivity  readings  and 
whenever  convenient  we  shall  replace  m  by  OK  )  and  vice  versa,, 


e.g, 


^rni  =  %  j  k 

In  (1):  Cfml  is  an  m  x  1  matrix.  Matrix  elements  ^  represent 
the  total  potential  for  a  combination  of  j  and  k  (which  is 


denoted  by  m) . 


is  an  m  x  1  matrix. 


Matrix  element 


is  given  by 


C 

2.7Tq£ 

and  may  be  interpreted  as  the  primary  potential  due  to  a  unit 
source  at  k  and  receiver  at  j  position. 

[  dmi"]  is  an  m  x  n  matrix.  Gmi  is  the  approximated 
geometric  factor  associated  with  cube  i  and  source-receiver 


(  ....  . 

--I 

•  -  J 


•  r  (  ) 

•  ;  •'•••  •  (  ) 
.  -  :  .  ...  ‘  ..  j 


"')  .  O.  . 


. 


# 


( 


) 


.  o  .■ 


.  .) 


■  'j 


i 


•  (  " .  o' 


combination  designated  by  m.  For  its  explicit  form  refer  to 
Appendix  I. 


CX.-l 


is  an  n  X  I  matrix. 


C 


Cf i‘  —  CTo  .  /* 

o'-  +  2  O'  > 


3-6  (2) 


represents  the  conductivity  contrast  factor  for  1’tt  cube. 
Approximation  (l)  assumes  that 


(a)  There  is  a  very  thin  layer  of  uniform  conductivity 

on  the  surface  of  the  earth.  Thus  O'j  -  °K  * 
a  constant,  where  Oj’  =  conductivity  at  receiver 
position  j ;  0^  =  conductivity  at  source  position 

k. 

(b)  The  electric  field  at  a  conductivity  interface  is 
just  the  primary  electric  field. 

(c)  The  form  of  Xias  gi-ven  by  (2)  does  not  follow 
from  the  approximation  itself.  It  was  suggested  by 
Vozoff,  and  C=  3.6  was  empirically  determined  by 
him  [  7].  It  may  be  mentioned  that  the  actual  Xt 
which  appears  in  the  approximation  is  given  by 


However,  the  behaviour  of  In  -gyj"  as  a  con¬ 
ductivity  contrast  factor  is  contrary  to  experience. 
It  appears  that  Xi  may  be  an  algebraic  function. 
Hence  the  empirical  suggestion  by  Vozoff  £7^ 

The  interpretation  scheme  (l)  is  based  on  the  fundamental 


hypothesis  mentioned  earlier.  The  hypothesis  is  moderately 
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successful  in  actual  practice.  The  reason  for  this  may  be 
seen  by  examining  the  physical  situation.  Suppose  we  have 
only  one  disturbing  cube.  According  to  the  hypothesis  we 
have  to  define  a  region  of  interest,  and  divide  it  Into  a 
regular  array  of  cubes.  In  the  solution  of  (l),  the  disturb¬ 
ing  cube  shows  up  as  such,  and  the  remaining  cubes  show  up 
as  either  having  conductivity  equal  to  or  very  slightly 

different  from  it.  This  behaviour  makes  it  possible  to  get 
a  reasonable  fit.  Since  (l)  is  not  exact,  a  situation  some¬ 
what  like  the  above  is  really  intuitively  suggested.  If  (l) 
was  exact  and  the  physical  situation  perfect,  then  and  only 
then  would  we  expect  the  disturbing  cube  to  be  the  only  one 
to  show  up  in  the  solution. 

In  applying  this  hypothesis,  several  practical  problems 
must  be  considered.  One  concerns  the  optimal  size  of  the 
blocks  into  which  the  subsurface  is  to  be  subdivided.  This 
should  be  governed  by  electrode  spacing,  and  should  vary  with 
distance  (horizontal  as  well  as  vertical)  from  the  electrode 
traverse.  Another  concerns  the  amount  of  information  avail¬ 
able  from  a  given  number  of  readings.  Thus  it  may  be  asked, 
that  for  a  given  number  m  of  readings  on  the  surface  of  the 
earth,  what  should  be  the  number  n  of  subsurface  blocks  of 
which  k  ^  n  are  anomalous,  so  that  all  the  anomalous  blocks 
may  be  completely  determined  and  vice  versa.  These  and 
similar  considerations  are  very  important,  and  obviously  will 
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have  a  great  bearing  on  devising  optimum  prospecting  plans. 

Vozof f  [6,7]  and  Ness  ^4  J,  have  discussed  some  of  these 
questions  in  detail. 

It  is  interesting,,  however,  to  look  at  the  last  question 
from  the  point  of  view  of  information  theory. 

Suppose  there  are  n  identical  coins  and  exactly  one  of 
these  is  false,  being  either  lighter  or  heavier  than  the 
good  coins.  It  Is  required  to  detect  this  coin  in  a  minimum 
number  of  weighings  m,  by  an  equal  arm  balance. 

In  reference  £8]  there  appears  a  solution  to  this  problem. 
The  solution  is  based  on  the  following  assumptions : 

1.  The  permutations  of  the  states  of  n  coins  are 
distinguishable.  This  can  be  done  by  labelling  the 
coins  as  1,  2,  . ...,  n.  (This  may  appear  to  conflict 
with  the  word  "identical’'  in  the  problem,  but  it 
need  not  bother  us  since  a  labelling  scheme  will  be 
used  in  the  application). 

2.  Any  one  of  the  coins  could  be  false. 

The  following  solution  is  given: 

m 

n  <  I 

z 

Thus  to  detect  a  false  coin  from  a  set  of  12  coins,  3  or 
more  weighings  are  necessary. 

Assuming  any  k  ^  n  -  1  coins  are  false,  the  above  formula 
is  extended  to  take  the  form 


J2L 

K 


z 
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Suppose  now  wo  have  a  group  of  n  cubes  embedded.  In  a  uni¬ 
form  medium  of  conductivity  ,  and  that  the  conductivities 
of  k  of  these  cubes  differ  from  .  The  conductivities  are 
designated  by  ,  2=0 *  lj  2,  ...^n,  and  we  can  use 

equation  (3)  to  relate  n,  the  number  of  cubes,  m  the  number 
of  observations  and  k  the  number  of  "anomalous"  cubes.  In 
the  resistivity  problem  usually  k  =  n,  and  in  all  we  have 
n  +  1  conductivities  at  ,  Z  =  0,  1,  ...  n.  Hence  in  (3) 
we  put  k  =  nj  and  replace  n  by  n  +  1  to  get 

ro 

ft  ^  l — (3a) 

Z 

(where  n  and  m  are  integers)  and  assume  that  it  gives  an 
indication  of  the  minimum  number  of  readings  m,  required  to 
determine  the  conductivities  of  n  regions. 

In  typical  situations  where  perhaps  30  to  50  observations 
have  been  taken  along  a  traverse,  one  should  not  expect  to 
be  able  to  resolve  more  than  about  one -third  of  this  number 
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III.  Errors  and  Their  Estimates 


1.0 


A  measured  number  by  itself  signifies  nothing  that  could 
be  safely  interpreted  theoretically.  If  an  aggregate  is  avail¬ 
able,  only  then  can  the  theoretical  significance  of  the  measure¬ 
ments  be  assessed.  For  in  that  case  only  have  we  the  means  for 
determining  the  error,  or  measure  of  precision,  of  results. 

A  "true"  value  of  a  measured  quantity  cannot  be  determined. 
What  we  accept  empirically  as  a  final  value  is  really  a  maxi¬ 
mum  likelihood  value.  The  concept  of  maximum  likelihood  attains 
meaning  only  if  a  sufficient  number  of  measurements  is  avail¬ 
able.  From  this  it  is  possible  to  plot  a  distribution  curve 
that  represents  the  Information  contained  therein.  The  proced¬ 
ure  more  commonly  followed  Is  to  assume  a  Gaussian  distribution 
and  to  apply  least  squares  fitting  to  determine  a  most  likely 
value.  The  main  reason  for  the  wide  usage  of  the  least  squares, 
method  is  its  simplicity  and  the  ease  with  which  It  can  be  ap¬ 
plied.  However,  it  must  be  pointed  out  explicitly  that  one 
cannot  always  apply  the  method  of  least  squares,  and  for  certain 
types  of  observations  the  method  will  not  determine  a  most 
probable  value,  and  indeed  may  determine  a  least  probable  value. 
An  example  illustrating  this  point  may  be  found  in  [5  Jon  pages 
217  -  218. 

In  resistivity  work  magnitude  determination  involves 
reading  the  coincidence  of  a  pointer  on  a  scale.  The  limita¬ 
tions  of  an  instrument  and  our  senses  necessarily  Imply  an 
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inherent  uncertainty  in  the  determined  magnitude.  These 
errors  may  be  termed  the  resultant  experimental  error. 

These  errors  involve,  really,  a  rounding  off  of  the  magnitude, 
and  are  also  termed  roundoff  errors. 

The  observed  data  are  usually  approximated  by  theoretical 
expressions  based  on  certain  hypotheses.  This  introduces, 
what  we  may  call  theoretical  or  truncation  error.  A  trunca¬ 
tion  error  corresponds  to  the  fact  that,  whereas  an  exact 
result  would  be  afforded  by  an  infinite  sequence  of  steps, 
the  process  is  "truncated"  after  a  certain  finite  number  of 
steps.  (indeed  the  linear  approximation  given  by  equations 
(l)  is  got  by  truncating  the  process  of  successive  approxima¬ 
tions,  after  a  first  approximation  has  been  made.) 

In  using  a  theoretical  approximation  to  extract  useful 
information  from  observed  data,  the  question  of  estimating 
errors  becomes  important.  The  useful  information  which  a 
theoretical  approximation  is  designed  to  extract  usually 
consists  of  numerical  values  of  certain  parameters.  In 
particular  these  are  the  numerical  values  of  X  i*  for 
7  =.  )j  2j  •  •  y  H  in  our  case.  The  question  arises.  If  a 

number  of  determinations  have  given 

<3,  ,  Cu,  •  •  ■  ,  Q-m 

as  values  of  a  certain  quantity  in  various  situations,  what 
can  we  say  (in  a  maximum  likelihood  sense)  about  the  parameters 


which  appear  in  the  theoretical  approximation  of  this  quantity? 

In  particular  if  were  the  observed  potentials 

then  what  can  we  say  about  ? 

It  turns  out  that  if  Q*n  are  given  without 

any  indication  of  the  way  in  which  they  were  arrived  at,  a 
definite  estimate  in  a  maximum  likelihood  sense  can  only  be 
made  under  certain  arbitrary  assumptions.  Prom  these  assump¬ 
tions  it  is  possible  to  find  a  law  of  distribution  of  errors. 

If  the  law  of  distribution  of  error  is  known  or  assumed, 
and  the  resulting  fit  is  not  good,  then  necessarily  we  must 
re-examine  the  theoretical  approximation.  It  is  however 
usual  that  the  law  of  error  is  unknown,  and  one  has 
to  assume  its  form. 

In  practice,  one  assumes  the  law  of  error  (usually  Gaussian 
Law  of  error)  then  develops  a  theoretical  approximation,  so  as 
to  arrive  at  a  good  fit.  Clearly  this  is  an  apriori  procedure, 
and  depends  largely  in  one*s  belief  in  the  assumed  law  of  error. 

In  the  following  some  error  estimates,  the  underlying 
assumptions,  and  the  corresponding  laws  of  error  are  stated. 

Also  we  shall  mean  by  "error"  the  difference  between  observed 
and  predicted  values  of  the  quantity  measured^ 

1.  If  we  assume  that  for  a  given  set  of  data,  the  best 

possible  estimate  of  the  measured  quantity  is  given  by 
minimizing  the  sum  of  the  squares  of  the  errors,  then  this 
assumption  implies  that  the  probability  density  of  the  errors 

*The  more  common  term  for  this  difference,  in  statistical 
analysis,  is  the  "residual." 
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is  given  by  an  expression  of  the  form 


where  y  is  the  probability  density  of  the  distribution, 
and  is  such  that  the  probability  of  an  error  in  the 
determination  lying  between  x  and  x  +  dx  is  equal  to 
y  dx.  h  is  a  constant  called  the  precision  constant. 
Fisher  has  shown  that  if  the  above  assumption  is 
strictly  obeyed  then  minimizing  the  sum  of  the  squares 
of  the  errors  gives  maximum  information  possible,  for 
a  given  number  of  readings.  He  has  further  shown 
that  if  the  above  law  is  valid,  but  the  sum  of  the 
absolute  value  of  errors  is  minimized,  then  approximate¬ 
ly  14$  more  readings  are  required  to  obtain  the  same 
accuracy  which  may  be  obtained  by  minimizing  the  sum 
of  the  squares  of  the  errors. 

Burnside  [il  has  shown  that  the  above  assumption 
is  reasonable  under  the  following  conditions: 

(i)  If  actual  error  is  made  up  of  a  sufficiently 

large  number  of  minor  errors  of  the  same  magni¬ 
tude,  each  of  which  is  equally  likely  to  be 
positive  or  negative,  then  it  follows  Gauss Js  Law. 

(ii)  Error,  if  it  arises  as  the  sum  of  a  sufficiently 
large  number  of  minor  errors  of  which  the 
probable  value  is  zero,  follows  Gauss *s  Law. 

(iii)  If  a  number  of  independent  component  errors 
all  follow  Gauss *s  Law,  and  if  the  resultant 
error  is  compounded  linearly  from  them  in 
any  way,  then  it  also  follows  Gauss* s  Law. 


It  is  possible,  however,  that  each  component  of 


resultant  error  may  not  follow  Gauss *s  Law.  Though 
this  law  is  very  popularly  used,  yet  one  must  exercise 
care  in  its  use. 

Assume  that  minimizing  the  absolute  sum  of  errors 
gives  the  best  possible  fit. 

Suppose  n  measurements  are  made  to  estimate  a 
certain  magnitude  whose  true  value  is  )(  .  Let 
these  n  values  be  a^,  ag  ...  an  arranged  in  ascending 
order  of  magnitude.  Assume  further  that  t  lies 
between  ar  and  ar+1.  The  sum  of  absolute  values  of 
error  I 


£  Is  least  for  \  =  ar 

£  is  least  for  )(  =  ar+-^ 

£  is  independent  of  ){ 


when  2r  y  n 
2r  n 


2r  =  n 


When  n  is  odd,  say  equal  to  2mr  +1,  £  is  least  when 


0  =  am?.  When  n  is  even,  say  2m*,  £  is  constant 

for  all  values  of  !(  from  am*  to  a^+1  and  this  value 
is  least. 


If  the  above  assumption  is  justified  then  the  law 


of  error  is  of  the  form 


c  . 
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and  it  can  be  shown  (Fisher:  2)  that  the  above 
estimate  is  such  that  for  a  given  number  of  readings 
maximum  information  is  obtained  by  minimizing  the 
absolute  sum  of  errors. 

Suppose  errors  in  a  system  of  equations  are 
denoted  by  ,  index  Z  referring  to  2^  equation. 

can  be  either  positive,  zero,  or  negative. 

We  can  express  €.z  as 


all  °iz  and  are  either  positive  numbers  or 
and  only  one  of  the  pair  can  non~ 


where 


zero. 


zero. 


The  above  shows  that  when  dealing  with  errors 
as  such,  we  can  have  an  estimator  of  error,  viz : 


Minimizing  ( (?  |  means  that  the  sum  of  absolute 


value  of  errors  is  minimized. 
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3.  Another  estimate  may  be  that  the  best  value  Is  got 
by  minimizing  the  2rvi"ti)  power  of  errors  Ui  1  • 


2 


i- 1 


2m-  l 

-  2f)  =  o 


It  can  be  shown  tn  ,  that  when  nn  >  I 


where  0 {  |  =  minimum  Ui 

i 

Qn  =  maximum  ^x 

t 

This  shows  that  this  estimate  makes  the  determination 
of  true  value  depend  more  on  larger  and  smaller 


values  observed 


1  .  '.‘i  O'  '  :  '.C  )  _  "  'j_' 

i  -  t.  ’j.3  . 


' 
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Each  datum  of  apparent  resistivity  is  got  by 
reading  voltmeters  and  ammeters,  used  in  the  instru¬ 
mentation.  Assuming  these  meters  to  be  properly 
calibrated,  we  may  expect  the  errors  resulting  from 
the  use  of  these  meters  to  be  Gaussian.  In  the 
calculation  of  the  geometrical  coefficients  G^  an 
error  of  the  order  of  1%  is  indicated.  Whether  or 
not  this  is  Gaussian  is  open  to  question.  Further 
it  is  also  not  quite  clear  if  the  error  introduced 
by  the  first  approximation  used  to  arrive  at  equations 
(l)  is  Gaussian.  This  means  that  the  least  squares 
fitting  of  the  data  may  or  may  not  be  the  best 
approach.  In  view  of  the  semi quantitative  results 
got  by  using  the  method  of  least  squares,  in  the 
following  an  attempt  is  made  to  use  a  fitting 
approach  in  which  the  sum  of  the  absolute  values  of 
the  errors  is  minimized.  Before  proceeding  any 
further  an  example  is  given,  showing  that  minimizing 
the  sum  of  the  squares  of  the  errors  may  not  always 
minimize  the  total  error  in  the  solution.  The 
total  error  is  defined  as  the  sum  of  the  absolute 
value  of  the  errors. 

In  the  example  the  data  -f(*)  were  required  to 
be  fitted  to  a  straight  line.  The  two  solutions 


are  shown  in  Table  I 
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TABLE  I 


X 

0 

1 

2 

3 

4 

Total  absolute 
error 

Given  qx) 

1.00 

3.85 

6.50 

9.35 

12.05 

L.S. 

1.03 

3-79 

6.55 

9.31 

12o07 

.20 

L.P. 

1.00 

3.76 

6.53 

9.29 

12.05 

00 

1 — I 

• 

L.S.  E 

Least  squares 

approximat i on 

to  /(*) 

L.P.  =  Linear  Programming  approximation  to  fj£) where 
sum  of  the  absolute  value  of  the  errors  was 
minimized. 
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IV.  Linear  Programming 

Linear  Programming  describes  the  solution  of  the  follow¬ 
ing  problem: 

Find  a  column  vector  X  =  (X-^  X2,  ...  Xn),  which  minimizes 

f  =  cX 

Subject  to  the  conditions 

GX  ^  b 

X  0  1 

where 

G  =  (G  )  is  m  x  n  matrix 

b  =  (bj  . . .  bm)  is  a  column  vector 

c  =  (c]_  .  .  .  cn)  is  a  row  vector 

This  is  termed  the  1 primal*  problem. 

A  solution  satisfying  I  is  termed  a  feasible  solution. 

A  basic  feasible  solution  (b.f.s.)  is  a  feasible  solution  with 
no  more  than  m  positive  Xj  .  A  b.f.s.  which  maximizes  f  is 
termed  the  optimal  solution. 

Dual  Problem 

Associated  with  every  primal  problem  is  the  *dual*  problem 
Find  a  row  vector  W  =  (w^  . „ .  Wm)  which  maximizes 


g  =  Wb 
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subject  to  the  conditions 

WCr  ^  C 

w  ^0  II 

Solution  of  the  primal  problem  provides  useful  informa¬ 
tion  about  the  solution  of  the  dual  problem.  In  fact  if 
X°  is  an  optimal  solution  to  I  and  W°  optimal  solution  to  II, 
then 

max  f  =  cX°  =  min  g  =  W°b. 

Standard  methods  exist  for  the  solution  of  L.P.  problems 
[9b  Perhaps  most  well  known  is  the  Simplex  Method.  However 
the  Revised  Simplex  Method  [31  is  more  suitable  for  computer 
calculations  and  will  be  used  for  the  numerical  computations 
in  this  work. 

The  success  of  L.P.  method,  especially  for  large  problems 
depends  on  the  use  of  fast  electronic  computers.  The  size  of 
an  L.P.  problem  which  can  be  solved  on  an  electronic  computer 
is  usually  limited  by  the  number  of  rows  in  G  matrix.  A 
problem,  when  formulated  as  a  primal  problem,  may  be  too 
large  for  the  computer.  However,  its  dual  problem  may  be  of 
acceptable  dimensions.  As  was  mentioned  earlier,  one  can 
extract  the  optimal  solution  of  the  primal  problem,  if  the 
final  iteration  in  the  solution  of  the  dual  problem  is  known. 
It  is  to  stress  this  point  that  the  two  problems  are  stated 
together. 

The  method  of  solution  of  L.P.  problems  is  iterative  in 
nature,  and  is  easily  adaptable  for  computer  programming  and 
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solution.  This  allows  one  to  follow  through  the  successive 
stages  of  the  solution. 


V.  Formulation  of  the  Resistivity  Interpretation 

Problem  as  a  Linear  Program 


The  resistivity  approximation,  for  unit  current,  may 
be  written  as  follows : 

[%]  -  [Q  - 

Putting 

^  [ft  -ft]  =  [  Ll 

we  can  write 

where  .  .  (3) 

pj-  -  <Jo 

+  zd0 

We  introduce  in  the  approximation  a  column  vector 
and  interpret  it  as  an  error  column  vector.  Thus  the 
uncertainty  in  is  removed,  and  we  can  write  a  system 

of  equalities,  instead  of  the  approximate  equations. 

has  m  elements. 
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I’he  elements  oi  ["x^Jand  can  be  either*  negative  or 
positive,  while  a  linear  program  and  Its  method  of  solution 
requires  that  we  deal  with  positive  variables  only.  We  can, 
however,  reduce  the  above  system  of  equations  to  the  canoni¬ 
cal  form  of  linear  program  constraints  by  expressing  the 
elements  of[Xj_]and  as  a  difference  of  two  positive 
numbers,  thus 

[x,i  --  [x;-xn,  [£„i  =  K-f.] 

Finally  we  could  write 

[G-wi] [ y,  -  X,-  ]  +  [°g~i  - 

[yH  [  x/] 

[  Z'  0  ^  ^,1  0 


and  minimize 


where  M  is  a  row  vector,  all  the  elements  of  which  are 
unity.  C<W]  has  m  elements. 

Ill  is  a  standard  linear  program. 


In  III  it  is  required  to  minimize : 

/ 


■f  -  +  ~  £  Gw  («J*»  -  £.  +  ft,) 


which  is  the  sum  of  absolute  errors.  In  doing  this  we  do  not 
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have  any  way  to  ensure  that  the  sum  or  total  positive  error 
Is  the  same  as  the  sum  of  total  negative  error.  However, 
we  may  want  to  design  our  procedure  to  equalize  these  errors. 
This  can  be  done,  using  L. P.  methods,  by  reducing  the  absolute 
numerical  value  of 


-  P„]  -  -f  (A) 

The  simplex  method  first  minimizes  MM  (if  not,  it 
can  easily  be  made  to  do  so)  say  to  a  value  fp.  Then  it 
starts  maximizing  ,  thus  reducing 

as  a  whole.  It  is  quite  possible  that  f  may  change  sign. 

The  method  of  solution  is  iterative,  and  we  can  easily  have 
the  computer  print  the  solutions  of  each  iteration.  Each 
value  of  f  printed  will  be  smaller  than  the  previous  one. 
Suppose  the  sequence  of  the  values  of  f  is  as 


f 


k 


f 


P 


fr,f 


pi 


where  fr  is  the  value  of  f  just  before  f  changes  sign  and 
f^  its  value  immediately  after  it  has  changed  sign.  The 
sequence  shown  Is  in  descending  order  of  magnitude.  Geometri¬ 
cally,  we  can  show  the  above  sequence  as 
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Then  optimum  f  Is 


We  stop  as  soon  as  f  has  changed  sign, 
given  by 


This,  however,  does  not  mean  that  the  sum  of  absolute 
value  of  errors,  which  is  given  by  [  finite  +  is  also 

minimized.  It  may  or  may  not  have  been  minimized.  By 
reducing  the  numerical  value  of  £  Cvi^f0^  -  we  have  only 

ensured  that  the  sum  of  errors  entering  the  solution  with  a 
positive  sign  (i.e.  sum  of  the  non-zero  elements  of  DM)  is 
as  closely  as  possible  equal  to  the  sum  of  errors  entering 
the  solution  with  a  negative  sign  (i.e.  sum  of  the  non-zero 
elements  of  M). 

Flexibility  of  the  methods  of  L.P.  is  easily  seen. 
Firstly,  in  the  case  when 

.f  =  CM  c  M 

there  is  always  the  possibility  of  properly  combining  the 
solutions  corresponding  to  fp  and  f^,  which  may  improve  the 
final  solution  considerably. 

In  cases  when  the  limitations  of  the  computer  do  not 
allow  a  complete  run  of  the  program,  we  still  have  quite  a 
few  ways  to  get  around  this  difficulty.  Firstly,  perhaps  a 
dual  formulation  of  the  problem  may  help.  Secondly,  the 
dual  formulation  failing,  we  may  run  the  program  in  steps. 
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For  convenience  of  writing  denote 
keeping  Its  significance  In  mind. 


Cxi-  xH 

by  M, 

Thus  we  may  1 

lirst  run: 

K-.e M  -M  =  [Kl 


and  then  run 


is  a  solution  to  the  complete  problem.  By  making  a  suitable 
choice  of  A  it  is  possible  to  reduce  the  sum  of  absolute 
value  of  errors  even  more. 

Fitting  of  data  by  successive  runs  is  quite  successful, 
and  by  devising  a  suitable  procedure  one  can  also  exercise 
control  to  some  extent,  on  the  solution  of  the  problem.  This 
is  quite  important,  especially  when  one  desires  to  judge  the 
physical  situation  subjectively  and  control  the  fitting  of 
data  accordingly.  This  is  possible,  because  it  is  a  feature 
of  L.P.  that  it  can  be  used  to  minimize  or  maximize  the 
effect  of  any  variable  or  variables. 

Minimizing  C^Vv'”!  L  ^  ^  R^') 

does  not  necessarily  mean  individual  errors  o(y.  and  (3^-  are 
also  being  minimized.  In  the  L.P.  process  of  reducing  the 
sum  it  is  possible  that  a  few  of  the  errors  actually  become 
larger.  However,  it  is  also  desirable  that  individual  errors 
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must  be  made  as  small  as  possible.  To  achieve  this  It  Is 
necessary  to  fit  the  data  In  a  systematically  controlled  way. 
For  example.,  It  was  found  that  a  problem  solved  by  linearly 
combining  the  solutions  of  the  two  problems 

r -rpj  - 

did  not  give  an  acceptable  fit  to  our  data.  The  reason  for 
discarding  the  solution  was  that  though  the  sum  of  the  absolute 
value  of  error  was  in  the  range  which  could  be  accepted, 
some  of  the  individual  errors  were  quite  large.  The  actual 
situation  was  that,  of  the  twenty-two  error  variables,  three 
or  four  variables  accounted  for  50  per  cent  of  the  total 
possible  error.  Accordingly  some  control  was  exercised  to 
get  a  ^smoother,,  fit.  This  was  possible  mainly  because  the 
fitting  was  done  in  successive  stages,  due  to  limitations  of 
the  computer. 

The  point,  however,  to  emphasize  is  that  it  is  quite 
possible  to  use  L.P.  methods  effectively.  Flexibility  of  LCP. 
methods  are  a  decided  advantage  when  it  is  desired  to  incor¬ 
porate  some  subjective  judgement  into  the  solution. 

In  formulating  the  problem  as  given  by  III  we  have 
introduced  (2n  +  2m)  variables,  and  we  have  only  m  equations. 
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and  (2n  +  2m)  inequalities.  The  variables  have  been  introduced 
as  follows : 

(l)  m  variables  in  .  These  are  those  error 

variables  which  enter  the  solution  with  a  positive 
sign. 


(ii) 


(iii) 


m  variables  in[j3m~^  .  These  are  those  error 
variables  which  enter  the  solution  with  a  negative 
sign. 

2n  variables  in  fxt  _  x:l 


In  any  solution,  however,  m  elements  of  [°4nJ  and  r 
will  be  zero,  since  of  °<K  and  only  one  can  have  a 

non-zero  value,  the  other  will  be  zero.  The  same  applies 
to  Xj,  X^,  only  one  of  the  two  corresponding  elements  will 
appear  in  the  solution,  the  other  will  be  zero.  Thus  the 
effective  solution  is  reduced  to  choosing  m  variables  out 
of  n  +  m  variables o  Since  n  will  be  greater  than  one,  the 
number  of  error  variables  in  the  solution  will  be  m  -  n. 

If  the  solution  progresses,  so  that  all  the  n  variables  con¬ 
tribute  to  the  solution  then  presumably  the  total  error  will 
be  the  least.  However,  as  also  mentioned  in  the  end  of 
Section  II,  it  will  be  necessary  to  limit  the  maximum  number 
of  cubes,  and  hence  the  number  n.  This  limitation  as  shown 
earlier  may  be  of  the  order  given  by  the  relation 
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The  introduction  of  so  many  variables  is  largely  due  to 
the  desire  of  solving  the  problem  with  minimum  of  subjective 
manipulation,  especially  where  signs  of  the  variables  are 
concerned.  However,  in  the  case  where  n  is  not  very  large, 
a  close  examination  of  the  problem  will  usually  indicate 
correctly,  which  of  Xzor  Xl will  enter  the  solution.  In 
other  words,  when  the  number  of  cubes  into  which  the  subsurface 
region  has  been  divided  is  not  very  large,  then  it  is  possible 
to  decide  by  inspection  of  the  equation  whether  a  cube  is 
resistive  or  conductive  relative  to  the  background  conductivity. 
This  assumes,  of  course,  that  the  signs  of  the  G  coefficients 
are  correct.  This  procedure  reduces  the  number  of  variables 
and  contributes  to  a  great  saving  of  time.  For  example,  in 
the  case  of  a  single  cube, 

Gn  xi  +  -  Pi  =  n 

If  the  sign  of  Gn  is  correct,  then  it  is  a  physical  necessity 
the  sign  of  G±1  Xj  be  the  same  as  b±.  This  fixes  the  sign  of 
X-^  uniquely,  and  one  of  (X^,  X”)  can  be  discarded. 
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VI .  Model  Data  Analysis 

Model  data  due  to  single  disturbing  cubes  are  analysed. 

We  denote  the  resistivity  contrast  factor  of  i ^  cube  by 
X^(i  =  1,  ...  n).  (Since  we  will  not  have  occasion  to  use 
the  X  axis  explicitly,  there  need  be  no  confusion.) 

The  model  experiments  were  carried  out  in  a  wooden  tank 
by  Dr.  Vozoff.  The  conductivity  of  saline  water  was  properly 
adjusted  and  a  resistive  or  conductive  cube  lowered  on  a 
stand  in  the  tank.  The  cube  was  of  2"  side.  Pole -dipole 
measurements  were  carried  out  with  the  dipole  spacing  being 

2".  r r] 

In  the  actual  L.P.  analysis  of  observed  data,  the  disturb¬ 
ing  cube  X^  was  "surrounded"  by  five  other  cubes  shown  in 
Pig.  2  and  distinguished  as  X,-,,  X^.,  X^,  X^,  X^. 

The  standard  graphical  presentation  of  data  is  shown  in 
Pig.  3  r  7^  Data  were  available  for  all  the  points  marked  © 
as  well  as  the  points  numbered  from  1  to  22.  The  data  refer¬ 
ring  to  positions  ©  were  not  used  mainly  because  the  coef¬ 
ficients  Gj_jk  were  not  available.  The  data  at  the  other  22 
points  marked  as  1,  2,  ....  22  were  ultilized.  To  each  of 
these  datum  points  there  corresponds  a  constraint  equation. 

We  shall  use  index  m  =  1,  2,  ....  22  to  denote  these  points 

and  the  corresponding  equations.  Thus  in  the  L.P.  Problem  we  have 
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22  equations.  The  combination  of  J  and  k,  and  hence  the  m  to 
which  these  equations  refer  are  tabulated  with  the  results  of 
the  analyses  of  the  data.  Only  J  Is  given.  It  being  assumed 
that  the  other  end  of  the  dipole  is  on  location  (j  -  l).  A 
typical  pole-dlpole  configuration  Is  shown  in  Fig.  3*  where 
the  receiver  is  at  (2,1)  and  source  at  (4).  This  combination 
of  source  and  receiver  positions  Is  assigned  an  index  m  =  1. 

Thus  this  combination  in  the  tables  I  -  VI  appears  as  j  =  2, 
k  =  4,  m  =  1.  In  the  figures  only  the  m  index  is  shown,  under 
the  black  dot. 

Resistive  Cube 

Following  are  the  results  of  analysis  on  the  data  when 
the  model  tank  had  a  single  resistive  cube  as  the  disturbing 
element . 

Due  to  the  ultimate  size  of  L.P.  problem,  it  was  decided 
to  successively  approach  the  solution. 

From  the  nature  of  approximation  being  used  it  is 
evident  that  the  cubes  to  which  Xi,  X^,  X^,  refer  are  the  most 
important.  Hence  in  the  first  step  a  solution  is  approximated 
using  only  these  cubes  and  the  L.P.  equations: 

[PC  =tu 

Here  we  are  forcing  all  errors  to  be  negative,  and  are  denoting 
them  by  [>ml  Further  along  in  the  solution,  when  both  positive 
and  negative  errors  are  allowed,  they  will  be  denoted  by  IXml 
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The  numerical  results  of  this  step  are  given  In  Table  II. 
It  will  be  noted  that  all  errors  are  entering  as  -pm.  This 
effectively  means  that  X±  has  been  assigned  its  upper  limit. 
For  a  better  fit  X±  must  assume  a  lower  value.  This,  however, 
will  require  that  some  errors  enter  the  solution  with  positive 
sign.  This  may  be  seen  as  follows: 

Suppose  =  x|  +  X^  Xq  >  x'i 

The  2  fit  equation  is 


Therefore 
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When  X-^ - 5*  X-^  in  a  succeeding  iteration. 


then 
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->(Gnxi  -  Pi)  =  <v 
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z  can  be  positive  or  negative. 

It  is  easily  seen  that  numerical  value  of  errors  can  be 
decreased  by  properly  choosing  X^.  Two  different  approaches 
were  used  to  choose  X^.  These  are  labelled  as  A  and  B. 
Approach  A  was  used  in  the  analysis  of  resistive  cube  data 
now  under  consideration.  Approach  B  was  used  to  analyse  the 
data  from  a  conductive  cube.  However,  for  the  sake  of 
continuity  both  of  these  approaches  are  described  below. 

A.  In  Step  I  of  the  resistive  cube  analysis  /-£ , 

iii  || 

are  zero.  It  is  easy  to  find  an  upper  value  of  X^ 
such  that  when  - ^  x{,  errors  corresponding  to 
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Fa ft 1 al  Solution,  Resistive  Cube 


k 

J 

m 

bm 

Pm 

4 

2 

1 

5.56 

1.20 

4 

1 

2 

3.33 

3.87 

4 

0 

3 

7-50 

6.97 

4 

-1 

4 

0.555 

0.220 

3 

1 

5 

2.78 

6.53 

3 

0 

6 

14.4 

10.  5 

3 

-1 

7 

2.22 

0.80 

3 

-2 

8 

1.39 

0.77 

2 

1 

9 

8.33 

0.00 

2 

0 

10 

27.8 

21.9 

2 

-1 

n 

7.22 

2.36 

2 

-2 

12 

0.555 

1.05 

2 

-3 

13 

1.11 

1.04 

1 

0 

14 

75.0 

36.1 

1 

-1 

15 

25.0 

12.1 

1 

-2 

l6 

4.44 

1.80 

0 

-1 

17 

66. 7 

32.6 

0 

-2 

18 

13.9 

16.2 

0 

-3 

19 

6. 67 

3.84 

-1 

-2 

20 

25.0 

10.0 

-1 

-3 

21 

5. 56 

6. 84 

^310  -  111 


xi  = 

0881 

resistive 

X3  - 

.012 

conductive 

x4  = 

.012 

conductive 

Total 

error 

entering  with  negative  sign 

*177. 

Pole -Dipole  Measurements. 

Pole  Source  at  k 

Dipole  Receiver  at  (j„  j-l). 


34 


Pq  and  p  assume  positive  values,  and  errors 
corresponding  to  all  other  B  are  still  negative. 

Thus 

^1  =  ^11^1  “  Pj_  G  1  =  9)  22 

^-0  1  *  1  ...  8,  10,  ...  21 

Next  we  equate  the  effects  of  cubes  corresponding 
to  X2,X^,  X^  to  the  major  G^  values.  This  has  two 
advantages;  time  of  the  second  run  is  greatly 
reduced,  and  possibility  of  errors  entering  as 
positive  errors  is  increased.  A  final  solution  can 
be  easily  deduced.  Thus  the  second  step  is 

g±2x2  +  g±5x5  +  G±£X£  "  =  Gnxi  ~ 

for  i's  chosen  from  (l,  ...  ,  22) 

In  the  actual  solution  for  the  resistive  cube,  Xq  =  .2 
was  taken.  It  was  noted  that  total  error  changed 
from  ^  177  to  ^70.  The  second  step  was  run  and 
results  are  shown  in  Table  III.  The  error  column 
has  been  labelled  as  G j_  since  now  both  positive 
and  negative  errors  appear  in  the  solution. 

B.  In  this  approach  we  choose  an  Xi  such  that  when 

,  ' 

X-^ - ^  X-^,  approximately  a  half  of  the 

assume  positive  values.  We  incorporate  this  change 
of  sign  in  the  original  equations.  Thus  finally  we 
have  a  set  of  equations  in  which  approximately  half 
of  error  variables  are  positive  and  other  half 
negative,  and  we  re-run  the  whole  program.  The  error 
will  be  greatly  reduced.  Residual  error  can  further 
be  reduced  by  fitting  more  cubes.  This  approach  is 
used  in  the  analysis  of  the  conductive  cube  data. 


TABLE  III 
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k 

j 

m 

bm 

rji 

4 

2 

1 

5.  56 

0.71 

4 

1 

2 

3.33 

0.10 

4 

0 

3 

7.50 

-1.73 

4 

-1 

4 

0.555 

0.08 

3 

1 

5 

2.78 

-1.51 

3 

0 

6 

14.4 

-0.91 

3 

-1 

7 

2.22 

0.62 

3 

-2 

8 

1.39 

0.20 

2 

1 

9 

8.33 

-0.52 

2 

0 

10 

27.8 

-4.08 

2 

-1 

ll 

7.22 

1.35 

2 

-2 

12 

0.555 

0.320 

2 

-3 

13 

1.11 

-0.17 

1 

0 

14 

75.0 

-2. 86 

1 

-1 

15 

25.0 

-1.37 

1 

-2 

16 

4.44 

1.00 

0 

-1 

17 

66. 7 

0.00 

0 

-2 

18 

13.9 

-2.25 

0 

-3 

19 

6.67 

-0.29 

-1 

-2 

20 

25.0 

-0.35 

-1 

-3 

21 

5.56 

-2.11 

-2 

-3 

22 

8.33 

0.00 

Xp  =  .681 

x2  =  .295 

X3  =  o  012 

Total  absolute 

£  310 

resistive  Xij.  =  .012 

conductive  X^  =  .145 

conductive  X£  =  .145 

error  entering  with  both 

^20 

conductive 

conductive 

conductive 

positive  and  negative 

sign  20. 


Pole -Dipole  Measurements  -  Pole  Source  at  k 

-  Dipole  Receiver  at  (  j j-l). 


Conductive  Cube 
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Following  are  the  results  of  the  analysis  of  the  data, 
when  the  tank  has  a  conductive  (graphite)  cube. 

The  data  were  analysed  in  two  steps. 

Step  I. 

Using  X^,  X/|  the  following  L.P.  problem  was  solved. 

C ^mil  T  -^i "1  -  [Pml  =  fbrnl 

Step  II. 

Approach  B  was  used. 

Results  of  Step  I  follow  on  the  next  page.  Table  IV. 

In  Step  I  Xj  =  .924.  It  was  noted  that  when  we  put 

Xx  =  xj  +  X^  =  .724  +  .2  i.e.  X^  =  .2, 

the  errors 

P  ;  5=  1,  2,  4,  7,  8,  9,  15,  18,  22 

S' 

remain  negative,  while  errors  corresponding  to 

pq;  q  =  3,  5,  6,  10,  11,  12,  13,  14,  16,  17, 

19,  20,  21 

change  sign  and  become  positive.  This  change  was  incorporated 
and  the  program  re-run.  The  error  was  reduced  from  ^  89  to 
^  39.  Thus  with  only  cubes  Xj_,  X^>,  X^  the  error  is  less 
than  9 We  still  have  Xg,  X^,  X^,  which  we  could  use  to 
bring  the  error  down.  This  was,  however,  not  done  in  the 
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Partial  Solution,  Conductive  Cube 


k 

J 

m 

OVUUULbl  VL' 

km 

Pm 

4 

2 

1 

2.78 

3.47 

4 

1 

2 

7.78 

2.05 

4 

0 

3 

14.2 

1.89 

4 

-1 

4 

0. 866 

0. 266 

3 

1 

5 

13.9 

0.5 

3 

0 

6 

25.0 

3.1 

3 

-1 

7 

2.50 

1.25 

3 

-2 

8 

1.67 

0.53 

2 

1 

9 

1 — 1 

• 

1 — 1 

1 — 1 

8.2 

2 

0 

10 

50.0 

7  •  3 

2 

-1 

11 

8.67 

0.00 

2 

-2 

12 

1.67 

0.28 

2 

-3 

13 

1.94 

0.25 

1 

0 

14 

97.2 

14.2 

1 

-1 

15 

36.1 

11.3 

1 

-2 

16 

4.16 

0.16 

0 

-1 

17 

77.8 

19.8 

0 

-2 

18 

29.2 

7.9 

0 

-3 

19 

9.44 

1.03 

-1 

-2 

20 

22.2 

3.0 

-1 

-3 

21 

16.7 

2.4 

-2 

-3 

22 

2.78 

0.00 

-440 

^89 

\  - 

.924 

conductive 

R  = 

.032 

resistive 

x4  = 

.032 

resistive 

Total  error  entering  with  negative  sign  zcQ 9 

Pole-Dipole  Measurements  -  Pole  Source  at  k 

-  Dipole  Receiver  at  ( j ^  j-l). 
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present  work.  In  the  final  tabulated  results,  the  error 
column  Is  labelled  as  6m  ,  showing  both  positive  and  negative 


errors  entered  in  the  solution  (Table  V) . 

Approach  A  could  be  very  useful,  if  it  was  required  to 
Incorporate  information  from  other  sources  Into  the  analysis 
of  field  measurements.  A  feature  of  this  approach  is  that 
the  time  required  for  computing  the  second  run  is  much  less 
than  the  time  required  for  computing  the  initial  run.  In 
case  of  small  problems  this  can  come  in  very  handy. 

Approach  B  is  to  be  identified  with  reducing  the  absolute 


numerical  value  of 


It  more  or  less 


gives  symmetrical  fitting  of  the  data.  A  feature  of  this 
approach  is  that  we  do  not  subjectively  change  the  value  of 
Xt,  but  this  results  in  its  becoming  more  time  consuming 
than  approach  A.  This  may  not  be  justified  in  most  cases. 

The  computation  time  using  the  L.P.  approach  is  just 
about  the  same  as  the  time  required  when  the  method  of  least 
squares  is  used.  For  comparison,  a  least  squares  solution 
has  been  done  on  the  conductive  cube  using  the  same  three 
contributors  that  were  used  in  the  L.P.  solution.  This  is 
shown  in  Appendix  II. 
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TABLE  V 

Final  Solution,  Conduct;  l  ve  Cube* 
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k 

j 

m 

km 

£  rri 

4 

2 

1 

2.78 

-2.20 

4 

1 

2 

7.78 

0.00 

4 

0 

3 

14.2 

1.4 

4 

-1 

4 

0 . 866 

-0.023 

3 

1 

5 

13.9 

2.5 

o 

0 

0 

6 

25.0 

2.7 

3 

-1 

7 

2.50 

-0.43 

3 

-2 

8 

1.67 

-0.08 

2 

1 

9 

11.1 

-4.1 

2 

0 

10 

50.0 

4.5 

2 

-1 

ll 

8.67 

1.73 

2 

-2 

12 

1.67 

0.12 

2 

-3 

13 

1.94 

0.20 

1 

0 

14 

97.2 

8.4 

1 

-1 

15 

36.1 

-1.5 

1 

-2 

16 

4.17 

0.70 

0 

-1 

17 

77.8 

0.00 

0 

-2 

18 

29.2 

-0.2 

0 

-3 

19 

9.44 

1.09 

-1 

-2 

20 

22.2 

2.0 

-1 

-3 

21 

16.7 

3.5 

-2 

-3 

22 

2.78 

-1.41 

*Compare  with  least  squares 

2*440 

2*39 

solution  in  Appendix 

II. 

Xx  =  .734 

conductive 

X3  =  .025 

resistive 

X4  =  .025 

resistive 

Total  absolute 

value 

of  error  entering  with  both 

positive  and  negative 

sign  39 

Pole-Dipole  Measurements  -  Pole  Source  at  k 

-  Dipole  Receiver  at  ( j „  j-l). 
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VII.  Discussion  of  Actual  Computations 


The  geometric  factors  G^j^  are  surface  integrals  which 
appear  in  the  linear  approximation.  These  were  calculated 
numerically  [l  ~]  ,  and  are  accurate  to  about  1%.  The  factors 
take  on  both  positive  and  negative  signs.  As  a  result  of  the 
approximation  involved  in  their  determination  there  are 
regions  (j,k)  where  their  zeros  do  not  coincide  with  the 

zeros  of  the  measured  data.  The  computed  values  of  G. 

1  jk 

refer  to  pole-pole  configuration  of  the  source  and  receiver. 
However,  pole -dipole  data  were  analysed  and  the  G^.^  were 
used  to  calculate  the  coefficients  for  the  pole-dipole 
configuration.  Accordingly  probability  of  a  geometrical 
coefficient  having  a  wrong  sign  was  increased. 

Using  only  X1  in  the  analysis  it  was  seen  that  the 
error  could  not  be  reduced  to  the  order  of  5$  -  10 hence 
the  basic  hypothesis  of  dividing  the  region  into  cubes  had 
to  be  evoked  in  greater  degree.  Now  another  interesting 
observation  was  made.  When  the  program  was  run  with  all  six 
cubes  available,  only  one  or  two  showed  up  in  the  solution. 

For  effective  reduction  of  error  it  is  necessary  that  the 
cubes  other  than  X1  should  have  finite  conductivities,  slightly 
differing  from  the  background.  If  they  show  up  as  having  the 
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same  resistivity  as  the  background,  the  corresponding  X  will 
show  up  as  zero  in  the  solution,  and  the  error  will  not  be 
decreased.  There  seems  to  be  a  reason  to  explain  this.  As 
far  as  the  mathematical  structure  of  the  system  of  equations 
is  concerned,  a  wrong  sign  of  the  coefficient  of  X  in  any 
equation  will  affect  the  whole  solution,  and  will  either 
give  X  a  very  wrong  value  or  will  simply  make  it  zero. 
Obviously  the  nonappearance  of  certain  X's  in  the  solution 
is  due  to  some  anomaly  in  signs  of  the  coefficients  of  X. 

That  some  anomalies  are  present  is  quite  obvious  and  in 
some  measure  may  be  attributed  to  uncertainty  in  sign  of 
coefficients  G^.^.  However,  it  was  also  noted  that  anomaly 
in  sign  usually  showed  up  when  the  spacing  between  the  source 
and  receiver  was  large,  i.e.  when  the  larger  volumes  are 
being  explored.  It  appears  plausible  that  perhaps  some 
surface  and  boundary  effects  come  into  play.  Thus  the 
discrepancy  in  sign  may  be  tied  in  with  the  size  of  the  tank 
and  perhaps  in  a  bigger  tank  these  would  not  appear  at  all. 

On  the  basis  of  the  above  arguments,  signs  of  four  to  six 
coefficients  were  arbitrarily  corrected  in  the  complete 
analysis  of  a  set  of  data.  It  was  also  noted  that  if  conduc¬ 
tivities  X^  and  were  forced  to  be  equal,  and  if  X^  and  X^ 
were  also  forced  to  be  equal,  the  number  of  sign  discrepancies 


was  reduced. 
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The  method  Is  quite  successful  in  semiquantitati vely 
analysing  the  problem.  The  semi quantitative  nature  of  Inter¬ 
pretation  arises,  not  from  the  method  of  fitting  used,  but 
from  the  nature  of  contrast  factor  Xp *  and  the  linear  approxi¬ 
mation  itself.  We  notice  that 


C  X, 


£  (TV  ~  (^0 
e/z  +  2  60 


£  ^  3-4 


is  not  a  very  sensitive  parameter  as  is  remarked  in  L7  J 
The  parameter  Xj_  approaches  1  as  conductivity  increases,  and 
approaches  -  ^  as  resistivity  increases.  In  case  of  conductive 
cube  a  value  of  .73  is  noted,  while  in  case  of  resistive 
cube  a  value  of-068  is  noted,  for  this  parameter.  The 
absolute  value  .68  . 5  thus  does  not  fall  within  the  allow¬ 

able  range  of  Xj_ .  However,  it  must  be  remembered  that  the 

parameter  d  ~ — ~zr  is  entirely  empirical  and  some 

GY  f  2G0 

deviations  from  it  are  expected  as  may  be  clearly  seen  from 
the  graph  used  to  determine  £=3.6  [7J  • 

A  tremendous  improvement  over  least  square  fitting  is 
not  noted.  The  work  of  Ness  [4l  shows  that  the  least 
square  method  of  fitting  at  best  gives  results  in  the  range 
of  5%  -  10 %.  The  results  achieved  by  L.P.  methods  are  clearly 
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in  this  range*  Hence  the  latter  methods  are  at  least  as  good 
as  least  squares  and  due  to  their  flexibility  may  even  claim 
an  advantage  over  the  method  of  least  squares.  However,  for 
large  problems  L.P.  methods  tend  to  consume  more  and  more  time. 

The  aim  of  achieving  something  better  than  5$  error  in 
fitting  is,  from  the  point  of  view  of  the  actual  errors 
involved,  not  justified.  First  of  all  an  error  of  about  1% 
exists  in  the  computed  coefficients  Then  there  are 

experimental  errors  which  are  of  the  order  of  (2-3)$  at 
least.  Finally  an  unknown  amount  of  error  is  introduced  in 

the  linear  approximation  and  in  the  empirical  parameter 

6^'  '  &o 


C 


< ri  +■  26^0 


These  errors  on  the  whole  are  sufficient 


to  limit  the  solution  to  a  5$  error  fit. 

However,  the  purpose  of  this  work  was  to  see  whether 
or  not  L.P.  could  be  used.  That  it  can  be  used  is  demonstrated. 


For  example,  compare  result  of  Table  V  with  that  in 
Appendix  II,  Table  VI. 
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IX.  Conclusions 


The  results  obtained  show  that  the  L. P.  methods  can  be 
used  for  the  Interpretation  of  resistivity  data„  From  the 
point  of  view  of  the  computing  time  and  accuracy  involved, 
these  methods  appear  to  be  as  good  as  the  method  of  least 
squares,  at  least  on  these  two  examples.  The  L.P.  methods  are 
quite  flexible  and  provide  means  of  subjective  control  over 
the  solution.  In  cases  where  the  latter  is  necessary,  these 
methods  are  more  powerful  than  the  method  of  least  squares, 
and  can  be  gainfully  employed.  At  the  present  stage,  however, 
the  choice  between  the  two  methods  would  perhaps  be  dictated 
by  the  amount  of  data  to  be  interpreted,  by  the  type  of 
information  desired,  and  by  the  amount  of  time  the  analyst  has 
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The  problem  of  resistivity  interpretation  can  be  formula¬ 
ted  as  a  problem  in  Theory  of  Games.  It  may  be  mentioned 
here  that  Theory  of  Games  and  Linear  Programming  are  closely 
related  to  each  other. 

We  may  consider,,  in  principle,  that  a  list  has  been 
drawn  up  of  every  possible  state  of  the  matrix  of  cubes  into 
which  the  subsurface  region  has  been  subdivided.  Each  cube 
can  exist  in  one  of  three  states  (assumed  equally  likely)  i.e. 
it  can  be  more  conductive  than  background  (  O £  >  &0  )*  less 
conductive  than  background  (  G'z  ^  )  or  same  as  back¬ 

ground  (  ov  =  %  ).  Thus  correspondingly  the  matrix  of 

n 

cubes  can  exist  in  3  states.  In  the  terminology  of  Theory 
of  Games,  each  state  of  the  matrix  of  cubes  is  termed  as  a 
pure  strategy  of  nature.  The  geophysicist  makes  m  measure¬ 
ments  on  the  surface  of  the  earth.  He  then  attempts  to 
predict  the  state  of  the  matrix.  His  predictions  involve  the 
state  of  the  matrix  of  cubes,  and  also  the  experimental  error, 
which  can  be  either  positive,  negative,  or  zero.  Thus  we  may 
regard  the  geophysicist  as  having  3n+m  strategies.  The  aim 
is  to  make  predictions  in  an  optimal  way.  Games  of  this  type 
are  discussed  in  "Introduction  to  the  Theory  of  Games"  by 
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McKinsey  (McGraw-Hill).  However,  it  is  obvious  that  the  above 
formulation  will  have  to  be  suitably  modified  before  a 
practical  solution  can  be  devised. 

Another  interesting  suggestion  is  due  to  Vozoff.  When 
the  graph  used  to  determine  C  =  3.6  in  [jJ  is  closely 

*  —  6*  o 


examined,  it  is  clear  that  the  parameter 


G'z'  -hzG'o 


is 


either  - 1  or  /  for  most  parts  of  the  range  of  the  rati< 


iZUl 


/V 
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It  shows  variation  only  in  a  very  short  range  when 
.  Thus  if  empirically  it  is  hypothesized  that 
is  either  one  value  of  the  set  of  numbers 


I  J  then  under  suitable  assumptions,  a  formulation 
in  terms  of  game  theory  may  be  possible.  Statistical  games 
and  testing  of  hypotheses  are  discussed  in  "Theory  of  Games 
and  Statistical  Decisions"  by  Blackwell  and  Girshick  (John 
Wiley) .  In  many  cases  resistivity  methods  are  used  to  decide 
whether  a  subsurface  region  is  resistive  or  conductive  relative 
to  background  without  reference  to  actual  magnitude  of  oY 
In  such  cases  the  above  formulation  would  be  justified. 

A  realistic  linear  programming  model  requires  accurate 
and  reliable  numerical  values  of  the  coefficients 
Hence  it  would  be  useful  to  study  the  behaviour  of  solutions 
when  these  coefficients  are  allowed  to  vary. 
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The  Algebraic  Expression  for  the  Geometrical 
Coefficients  G^j^  and  Some  Numerical  Values 
Obtained  Therefrom 


In  the  resistivity  approximation  coefficients  G^j^.  for 
pole-pole  configuration  are  given  by 
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(See  Fig.  4) 


FIG. 
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GEOMETRICAL  PICTURE  FOR  A  SINGLE  CUBE 


The  coefficients  Gj_  were  computed  numerically  on  an 
electronic  computer  cn-  Prom  these,  the  geometrical 
coefficients  corresponding  to  pole-dipole  configuration  were 
computed  for  the  cases  when  the  centre  of  the  disturbing 
cube  was  1  unit  and  2  units  below  the  surface  of  measurement. 
These  numerical  values  are  shown  graphically  on  the  following 
page.  (See  Pigs.  5  and  6.) 
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Appendix  II 

Least  Squares  Solution,  Conductive  Cube 

Model  data  due  to  a  single  conductive  cube  (See 
Table  V)  were  analysed  by  the  method  of  least  squares.  The 
results  of  this  analysis  are  shown  in  Table  VI.  The  results 
obtained  by  the  two  methods  (Table  V,  Table  VI )  are  in  the 
same  range  of  accuracy. 
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Least  Squares  Solution,  Conductive'  Cube 


k 

J 

in 

km 

4 

2 

1 

2.78 

-2. 4l 

4 

1 

2 

7.78 

0.18 

4 

0 

3 

14.2 

-0.4 

4 

-1 

4 

0.867 

0.060 

0 

1 

5 

13.9 

4.6 

3 

0 

6 

25.0 

2.7 

O 

0 

-1 

7 

2.50 

-0.34 

3 

-2 

8 

1.67 

-0.13 

2 

1 

9 

11.1 

-3.0 

2 

0 

10 

50.0 

4.5 

2 

-1 

11 

8.67 

1.22 

2 

-2 

12 

1.67 

CM 

1 — 1 

0 

0 

1 

2 

-3 

13 

1.94 

0.01 

1 

0 

14 

97.2 

5*  6 

1 

-1 

15 

36.1 

-1.0 

1 

-2 

16 

4.17 

0.37 

0 

-1 

17 

77-8 

-3.0 

0 

-2 

18 

29.2 

1.10 

0 

-3 

19 

9.44 

0.81 

-1 

-2 

20 

22.2 

0.20 

-1 

-3 

21 

16.7 

3.74 

-2 

-3 

22 

2.78 

-1.92 

£  440 

£  36 

X^  =  0.746  conductive 

X^  =  0.020  resistive 

X4  =  0.020  resistive 


Total  absolute  value  of  error  entering  with  both 
positive  and  negative  sign  36 

Pole-Dipole  Measurements  -  Pole  Source  at  k 

-  Dipole  Receiver  at  (j,  j-l). 
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